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Abstract. Heisenberg's unsolved fundamental equation of the 
universe [3 US] has a coupling constant I which has the dimen- 
sion of length [L] . We consider a linearized version of Heisenberg's 
fundamental equation which also contains a coupling constant I 
with the dimension of a length and we solve this equation in the 
framework of a relativistic quantum field theory with a fundamen- 
tal length £ in the sense of our recently developed theory [5] and 
show that then one has £ = 1/{'\/2tt). This is done in two parts. 
In this first part we use path integral methods (and nonstandard 
analysis) to calculate all Schwinger- and all Wightman- functions 
of this model, as tempered ultrahyperfunctions and verify some of 
the defining conditions of a relativistic quantum field theory with 
a fundamental length, FLQFT for short. As an important inter- 
mediate step the convergence of the lattice approximations for a 
free scalar field and for a Dirac field is shown. 

The second part completes the verification of the defining con- 
ditions of FLQFT and offers an alternative way to calculate all 
Wightman functions of the theory. 
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1. Introduction 

1.1. Heisenberg's fundamental equation. The basic relativistic equa- 
tion of quantum mechanics called Dirac equation 

h d 

?-7m^ — V'la;) - rmplx) =0, Xq = ct, Xi = x,X2 = y,X3 = z (1.1) 

C OXfj^ 

contains a constants c (velocity of light) which is the fundamental con- 
stant in relativity theory, and Planck's constant h = 2Tih which is the 
fundamental constant in quantum mechanics. The dimension of c is 
[LT-i] and that of h is [ML^T"^]. W. Heisenberg thought that a fun- 
damental equation of Physics must also contain a constant I with the 
dimension of length [L]. If such a constant / is introduced, then the 
dimensions of any other quantity can be expressed in terms of combi- 
nations of the basic constants c, h and /, e.g., time [T] = [L]/[LT~-'^], 
or mass as [M] = [UVT-^] / {[UT-^]\L]) . 

In 1958, Heisenberg and Pauli introduced the equation 

h d 

-7/.^V^(a;) ± l^lt.l5Hx)Hxh^l5^i^) = 0, (1.2) 

which was later called the equation of the universe and studied in [51 [T^. 
The constant / has the dimension [L] and is called the fundamental 
length of the theory. 

Unfortunately, nobody has been able to solve this equation. At 
present, even in the more advanced framework of ultra-hyperfunction 
quantum field theory, we do not see how this equation could be solved. 
Accordingly we study a linearized version of this equations which in- 
herits the important property of a fundamental length / and which 
first has been studied by Okubo |22]. This linearized version is solv- 
able in the sense of classical field theory, i.e., the classical fields (plx) 
and ip{x) = ijj'{x)e'^^ '^^^^ solve this system when is a solution of the 
Klein-Gordon equation and ip' a free Dirac field of mass m. We write 
it in the form 

— j (l){x) = 



c dx'^ J 

and propose to solve the quantized version of these equations in the 
framework of a relativistic quantum field theory with a fundamental 
length as proposed recently by the authors [2] by constructing the 
Schwinger functions of the fields 0(x) and ip{x). And we do so by 
invoking nonstandard analysis and path integral methods. Thus we 
calculate the Schwinger functions by means of path integrals on the *- 
finite lattice with an infinitesimal spacing. As a result, the Wightman 
functions (i.e., the Wick rotated Schwinger functions) of the field ip{x) 
are not tempered distributions, but an tempered ultra-hyperfunction. 
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In the following we will work with the natural units c = h = 1. Then 
the system of equations f 1 1.3 1) reads 

(□ + m2)0(x) = O (1.4) 

{i-ff^d^ - m)ij = 2l^-f>'ij{x)(j){x)df,(j){x). (1.5) 

and they are the field equations of the following Lagrangian density: 

L{x) = Lpfix) + Lpbix) + Lj{x), (1.6) 

Lpfix) = ij{x){t^^d'' - m)iP{x), (1.7) 

Lp,{x) = ^{{d'^^{x)f-m'<P{xf}, (1.8) 

Lj{x) = 2l^{ij{xh^ij{x))(j){x)d''(l){x). (1.9) 



1.2. Relativistic quantum field theory with fundamental length 
(FLQFT). As indicated above we are going to show that the system 
(11. 4p - (11. 5p can be solved in the framework of a relativistic quantum 
field theory with a fundamental length (FLQFT) as developed in [2]. 
This theory is essentially a relativistic quantum field theory in the sense 
of Garding and Wightman [21] in terms of operator-valued tempered 
ultra-hyperfunctions instead of operator- valued tempered Schwartz dis- 
tributions. The localization properties (in co-ordinate space) of tem- 
pered ultra-hyperfunctions (for a technical explanation we have to refer 
to [2]) are very different from those of Fourier hyperfunctions and (tem- 
pered) Schwartz distributions. Tempered ultra-hyperfunctions distin- 
guish events in space-time only when their distance from each other 
is greater than a certain length i (A heuristic explanation of this 
property is given in [2]). In contrast to this, Fourier hyperfunctions 
and Schwartz distributions form a sheaf over space-time and thus ex- 
hibit essentially classical localization properties. On the other side the 
Fourier transforms of tempered ultra-hyperfunctions have essentially 
classical localization properties in energy-momentum space. Accord- 
ingly, compared with relativistic quantum field theory in the sense of 
Garding and Wightman (abbreviated as QFT), it is the locality condi- 
tion (condition of local commutativity) which needs a new formulation 
in FLQFT. Based on the notion of carrier of analytical functionals we 
proposed and used in [2] the notion of extended causality or extended 
local commutativity. 

With this notion of extended local commutativity a full set of defin- 
ing conditions for a relativistic quantum field theory with a fundamen- 
tal length has been given and such theories have been characterized 
in terms of a corresponding full set of conditions on their sequences 
of vacuum expectation values (n-point or Wightman functionals). In 
addition an explicit model for such a theory is constructed in [2]. This 
model is the (Wick) exponential of the square of a free massive field 0, 
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i.e., the field 

n=0 

The two-point functional of this field is formally 

in,pix)piy)n) = [l-Ag'Di~\x-y)r'^' 

where Dm \x — y) is the two-point functional of the field 0, and the 
fundamental length of this model is 

The major achievements of QFT are the proof of the PCT theorem, 
the relation between spin and statistics and the existence of a scattering 
matrix. In FLQFT the PCT and the spin-statistics theorems and the 
existence of a scattering matrix have been proven too. 

1.3. Motivation for FLQFT. Very briefiy we recall our motivation 
for our version of a relativistic quantum field theory with a fundamental 
length. 

The first question one has to answer is on which level of the theory 
the fundamental length should be realized. 

The established answer to this question is that the fundamental 
length should be realized on the level of the geometry of the underly- 
ing realization of space-time and accordingly the 'standard' approach 
to a (quantum) field theory with a fundamental length is to invoke 
non- commutative geometry [U [261 E] • 

We think that it is important to keep as many of the established 
physical concepts and results based on the traditional realization of 
space-time as possible and accordingly have proposed in [2] to realize 
the fundamental length on the level of the primary dynamical quanti- 
ties of the theory, namely the fields. In this way we can rely directly 
on the established physical principles (of field theory, relativistic co- 
variance, physical energy- momentum spectrum, quantum physics). As 
pointed out above then the only change necessary is that of the re- 
alization of the locality principle of standard QFT (when the type of 
generalized functions to be used in this theory is set to be tempered 
ultra-hyperfunctions). In this way we arrive at a relativistic quantum 
field theory in which the fundamental length is realized through special 
localization properties of the fields and in which the major achieve- 
ments of standard QFT are still valid. 

In the second part where we actually prove these localization prop- 
erties for our solution we give a brief technical explanation of the lo- 
calization properties of tempered ultra-hyperfunctions (see subsection 
1.2). 



2. Path integral quantization 

As announced we quantize this model by path integral methods. 
Formally, the time-ordered two point function is calculated as (see [H]) 



i>a{xi)i/Ji3{x2) expi |y Lj{x)dx | dV{ip,ip)dQ{(t)) 

= expi < / LFb{x)dx \ \ \ d(f){x) 

All these integrals have a rigorous meaning if the continuum space-time 
is replaced by a lattice. We will control the transition from the lattice 
to the continuum limit by methods from non-standard analysis. 

For positive integers M,N define L = MN and A = -y/vr/Af. Then 
the lattice F = r(M, A^) is 

T = {t = jA;j eZ, -L<j<L}. 

The lattice version of the differential operator —A + m'^ on = M^'^-^ 
is the following difference operator on F^: 

- A + : R^* 9 <l>(x) 

where is the vector of length A parallel to the /i-th coordinate axis. 
Let dG{^) be a Gaussian measure on W^'"^^ defined by 

(2.1) 

where C is the normalization constant such that J dG{^) = 1. Note 
that the exponent of this measure is the (Euclideanized; x° — > — 
X ^ y) discretization of the Lagrangian / Lph{x)dx. For later use 



we recall the following well-known formulae for Gaussian integrals on 
(see [H]). 



(27r)-"/Vdet A / e^^^'^) 



exp 



~{x,Ax) 



dx = exp 



(A) 



6 



S. NAGAMACHI AND E. BRUNING 



(2.)-»/^v/d^/,.,A.)exp 



dx = Tr(AA-^) (B) 



where Re A is strictly positive-definite and A is an arbitrary matrix. 
Note that the path integral on the finite lattice is the usual integral. 

Using (B), we calculate the covariance of the measure dG{^). We 
define a function 6{y) on by 6{y) = if y = otherwise 6{y) = 0, 
i.e., S{y) = A~^(5o,j/. Then (— + m^)6{x — y) is the kernel func- 
tion of the operator —A -|- and (— A^ + m^)5{x — y)A^A^ corre- 
sponds to the matrix A of the formulae (A) and (B) since the sum- 
mation is always accompanied by A^. The inverse matrix A^^ 

corresponds (— A^, + m?)~^5{x — y) (note that there are no additional 
A). In fact, (-eA^ + m^)6{x - ?/)AW for e = is m^6{x - y)A*A^ 
= m^5o.a;-j/A~"^A'*A^ and its inverse is m~'^6o,x-yA~^ = m~'^6{x — y). 
Now we can calculate the covariance of dG{^). 

^y^)^y2)dGi<^) = A^-^,^^ = (-A + mY^Vu 2/2) = SUvi - ^2). 

Using the lattice Fourier transformation, Sm{yi — 2/2) is representable 
as follows: 

-1 



3 

^(2-2 cos A) /A^ + 

.M=0 



5„(yi-i/2) = (2vr)-^^e*^(^-^^) 

(2.2) 

where the dual lattice F is given by 

T = {s = jri;j eZ,-L<j< L}, rj = V^/N. 

It converges, as M, N 00, (see the following section) to the two point 
Schwinger function 

SmiVi - y2) = (27r)-^ [ e'^(^-^^) [p' + m'] d^p. (2.3) 

of a neutral scalar field of mass m. 

In order to deal with the fermion field \1/ in the system fll.4|] - fll.Sp we 
need to do integration over Grassmann algebras, see pQ. Accordingly 
we define a measure dDl"^^, \1/^) on the Grassmann algebra generated 
by (see P) {^i{y), ^Uy); a = 1, . . . , 4, 2/ G F^: 

Ufld^ijMd^liy). (2.4) 

ygr* 0=1 
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where C is another normahzation constant, and 

vl/i = (vl/},...,vl/if, ^^ = {^1...,^IY. 
The matrices 7^ are related to the Pauh matrices (Tj by (j = 1, 2, 3) 

E f o-Q \ E ■ f ^ -io-j 

7o = 70 = n ' = -'^J ~ 



1 \ f l\ f -i\ f 1 



and the operators as discrete versions of the corresponding partial 
derivatives 9^ are defined as follows: 

V7 VI/ - / K'^'^iy) = ('^'^(y + ^m) - *fc(y))/A if ^ = 1, 2, 

^M^'^ \ V;vl/,(2/) = (M/,(y)-M/,(y-e^))/A iffc = 3,4; 



namely, 



= P+Vt + P-V- P± = (1 ± 7o'')/2. 



The idea to replace the partial derivatives in the continuum case by the 
forward-, respectively backward difference on the lattice as described 
above, has originally been developed in |21j . 



Remark 2.1. It is well known that the free fermion theory on the lattice 
defined by the action 

J2 *'(^) ( E7f [^'(a; + e,) - - e,)]/2A + m^\x)\ A^ 
xGr4 \ii=o J 

(2.5) 

suffers from the doubling problem. Wilson [25] has overcome this prob- 
lem by adding the term 

- J2 *'(^) ( + e/.) + ^'(a: - e^) - 2^i(x)]/2A J 

to (1231) 

It is also known that the doubling problem is due to the replace- 
ment of the partial derivative by the central difference (\l/(x -|- e^) — 
\l/(x — e^))/2A. If we replace by the forward difference respec- 
tively the backward difference V~ as we have suggested above, we have 
no doubling problems. Concretely, this is implemented in the Fermion 
Lagangian density (12.51) by choosing the forward difference for the com- 
ponents \l/2 of the Fermi field while the backward difference is used 
for the remaining components \l/3, \l/4. 

Kogut and Susskind [H] replaced the derivative of the space variable 
by half of the central difference, i.e., by (\l/(x-|-e^/2) — \l/(x — e^/2))/A. 
Then the doubling problem disappears but we must introduce the even 
lattice Fe and the odd lattice Fq and assign the subset (Fg U Fq) x 
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(Fg U To) X Fe) of (Fe U ToY to each field component as its domain of 
definition. For furtlier details about lattice fermion see 



In Section 4 we are going to show that the continuum limit of the 
covariance (two point function of the lattice Dirac field) 

3 

.M=0 



(2/1,2/2) = I ^l{yi)^}{y2)dD{^\<i/') (2.6) 

a,f3 

coincides with the Schwinger function Rm-a,p of the free Dirac field of 
mass fh: 



a,l3 



where 



Remark 2.2. Though nobody seems to doubt the convergence of the 
lattice approximations S^iyi ~y2) respectively 'R.m;a,i3{yi — 2/2) to their 
standard continuum forms Smiyi — 2/2) respectively Rm]a,i3{yi — 2/2) we 
could not find a proof. Maybe these convergence proofs are considered 
to be too tedious, especially in the case of Fermions due to the doubling 
problem. 

In Section 3, we prove this convergence by using nonstandard analy- 
sis [231 E], that is we show that for any infinitely large M, G *N, the 
standard part of Smiyi - 2/2) is Sm{yi - 2/2)- 

Remark 2.3. Sometimes, the proof by nonstandard analysis is sim- 
pler and clearer than the standard proof. For example, in order to 
prove lim^^oo f{n) = 00, for a function / : N ^ N, we must show: 
\fM3N\/n{n>N =^ f{n) > M). But in nonstandard analysis, 
we can use the formula I{x): x is an infinitely large number, and we 
have only to show Vn (/(n) ^ I{f{n))). The number of quantifiers 
is reduced in this nonstandard proof, and thus it is simpler and clearer 
(in technical terms: the syntactic complexity of the formula is reduced 
from a formula to a Hi formula. See |27j). 

In Section 4, the continuum limit of 'R.m;a,i3{yi — 2/2) is shown to 
be Rrh;a,i3- As anuouuced our prescription for avoiding the doubling 
problem works well here. 

Note that in sections 3 and 4, convergence of Schwinger functions 
is meant not in the sense generalized functions but in the sense of 
functions. 

Certainly, readers can skip sections 3 and 4 if they know or accept 
that these lattice approximations converge to their expected continuum 
limit. 
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Our overall strategy is as follows: 

(1) Construct the Schwinger functions 

[1 - 4Z^5^(yi - 2/2)T'^'^™(yi - y2) 

in the nonstandard universe; 

(2) by taking the standard part or continuous limit, we get stan- 
dard Schwinger functions [1 - Al'^Smivi - |/2)^]~"^^^-Rm(l/i - Z/2)- 
Sections 3 resp. 4 treat the continuous limit of Smiyi — ^2)^ 
resp. Ufniyi - 2/2). 

(3) by Wick rotation, we try to define Wightman functions 

lim [1 - AfDi~\xQ - le, x)^]^^/^ = lim [1 - 4/^5^(zxo + e, xf]~^/^ 

from Schwinger functions. But unfortunately, for this to give a 
mathematically well defined generalized functions e cannot be 
too small, actually e must be greater than i = //(v^tt); 

(4) in this way, the Wightman functions cannot be a tempered 
distribution but they can be ultra-hyperfunctions which satisfy, 
as we will prove later, axiom (RO), and i is the fundamental 
length according to axiom (R3). 

Next we describe our strategy of how to deal with the interaction 
in this model. We define the Euclideanized lattice Lagrangian density 
Li{y) which corresponds to the interaction Lagrangian Li{x) in (II. 9p 
as follows: 

-Ljiy) = ^'^iy)e^^'^^y^'j2^^ 

If we replace the differences in this definition by the corresponding par- 
tial derivatives (continuous limit) the above Lagrangian density Li{y) 
becomes the Euclideanization {x^ — » —iy'^, x ^ y) of iLj{x) as given 
in (O). 

Now we calculate the lattice version of the Schwinger functions of 
the interacting fields. The two point Schwinger function is 



*i(z/i)^J(y2)exp 5^L,(Z/)A4 rfD(vl/i,vl/2)rfG($) 

X ||exp (^^Lj{y)A'^ (iD(v^\ v^^)ciG($) | . (2.7) 
If we change the variables 
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then (12. 7p becomes 



m'\y{)^'\y2)dD{^>'\^>''' 



dG{^). 



As we are going to show the continuum hmit of 



is the two point Schwinger function Rm;a,i3{yi — 2/2) of the free Dirac 
field. Abbreviate h± = e^*'^/^/ and observe that the characteristic 
function of /i+$(?/2)) is 



Jth-<S>{yi)+ish+<S>{y2) 



dG{^) 



= exp --{th^Smiyi - yi)th^ + sh+Sm{y2 - yi)th_ 

+th^Sm{yi - y2)sh+ + sh+Smiy2 - y2)sh+} 

= exp -^{2tsl^SUyi - ^2) - ^t'/'5„(0) + isH^S^O)}. 
By using the relation (see formula (A)) 

we find 

= W J dtdse-'"l^e-'"'^ j e**V2^-*fei)e*v^^'^+*(^2)dG'($) 



= [(1 - 2z/25„(0))(l + 2z/25„,(0)) - 4/X(yi - ^2)'] 

where we used formula (A) of Gaussian integration for y = 0. 

The value of the two-point Schwinger function at the origin in the 
lattice approximation diverges in the continuum limit, i.e., 5^(0) = 
5m(0; iV, M) ^ 00 as iV, M ^ 00. In fact. 



^(2 -2cosp^A)/A2 + 

.M=0 



m 



m 
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-1 



-.(27r)-M ^4|p|V7r2 + m2 d'^p {N ^ oo) 

— > oo as (M — > oo). 

As usual we eliminate this divergence by interpreting the above prod- 
ucts respectively power series in the sense of Wick products. A compact 
form to define Wick products is as follows (see [TU]): 



e^^.^.^... {ith±^y)f : /n!] = e^'^ ^ — e 



t^fiSm{0)Jth±'i>{y) 



n=0 



Then we have 



and 



= (27r)-i j dtdse-'"/^e-'"/^ j : e'^'^^-'^^y^^ : : e'^'^+^^y^^ : dG{^) 
= (27r)-i J dtdse~''/^e-''/^ exp - {2ts/X,(yi - 1/2)} 

= [1 - Ai'^sUvi - y2)T'^' . 

Thus the two point Schwinger function of the field ip in lattice approx- 
imation is 

[1 - 4/^5^(1/1 - y2YY^^^ 'R'fh;aAyi - ^2), 

and its continuous limit (Section 4) is 



In order to construct the complete theory the system of Schwinger 
respectively Wightman functions of all orders n E N has to be con- 
structed. We show here how all ra-point functions of the interacting 
fields and ijj can be calculated in the lattice approximation. For the 
n-point Schwinger functions of the fields : g**j '?^(^) ; gn^j 



/n 



g ~tjt^,hrjhrf,Sm{yj-yk) 

j = l j<k 

where rj = + or = — , and 



/n 

» n „ n 
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i=l j<k 

where the matrix C = {cj^k) is given by 

Cjj = 1, Cj-fc = Ckj = 2hr^hrj'^Smiyj - Vk) 3 < k. (2.8) 

where again formula (A) has been used. 

Similarly, introduce the matrix A = (aj^k) by 

ttjj = 1, Cj-fc = akj = 2hr^hrj'^D\^\xj - Xfc), j < k, (2.9) 

Dl;;\xo,x) = Sm{iXo,x). (2.10) 

Then the ra-point Wightman function 

(0|p«(xi)---p('^)(xO|0) 

of the field 

p(^)(x,)=:e-(-i)''^^''<^(^^)': (2.11) 
is the Wick rotation of the Schwinger function (12.81) . i.e., 

(det A)"i/2_ 

This field has been studied in some detail in [12], see also 



Next, let ipo{x) be the free Dirac field of mass m and introduce the 
field components ip^{x) = ipoi^x), ip'^{x) = i^oix). Denote the Wight- 
man function of the free Dirac field ipo{x) by 

w^; jxi, . . . , X.) = {n, ra\ (x,) ■ ■ ■ vz i^rM 

and let S'q ^.(yi, . . . , ?/„,) be its Schwinger function, where r = (ri, . . . , r„), 
a = (ai, . . . , a„). Then the n-point Schwinger function of ipiy) is 

{det Cr'/'Sl,{yu...,yn), (2.12) 

where Sm{yj — yk) of (12. 2p is replaced by its continuous limit Smiyj — yu) 
(12. 3p and where the matric C is defined in 



Similarly, introduce the matrix A = (aj^k) as in (12.91) . Then the n- 
point Wightman function W^(xi, . . . , x„) of the field ip{x) is the Wick 
rotation of the Schwinger function (I2.12p . i.e., 

W;(xi,...,x„) = {detA)-'/^W',Jxi,...,Xn). (2.13) 

In Section 5, we show that the Wightman functions of ■iIj{x) are 
not tempered distributions but tempered ultrahyperfunctions which 
are studied in [H], [181 [2] . The axiom (RO) of reference [2], modified for 
the case of Dirac fields is verified. 

In part II of our investigations of this linearized model of Heisenberg's 
equation [3] , it is shown that the present model satisfies all the axioms 
of relativistic quantum field theory with a fundamental length. 
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3. Convergence of the lattice approximation of the two 
point functions for free scalar fields 

For positive integers M, put L = MN and let F be the 1- dimensional 
lattice 

F = {x = jA;jeZ, -L<j<L}, 
with spacing A = y^/M. Its dual lattice 

f = {p = jri;jeZ, -L<j<L} 

then has the spacing rj = y/ii/N . Let be the vector parallel to the 
/i-th coordinate axis with length A, and the forward respectively 
backward difference in direction defined by 

Then we have 

V+e'P" = e'^'" = zg^e**'" (3.1) 

V;e*P^ = e*P^ = ig^e^P^ (3.2) 

where = (1 - e'^P'''^) / (iA) . Note that 

S7^S7;e'P^ = -\q,\'e'P^ = _ ^ - 2 cosp^A ^,^,^ 

Accordingly we define a linear operator —A+m^ = — Yl'l=o ^^'^Ji'^^'^ 
on M'"^ = M^'^^ (second order difference operator on the lattice F^) by 

- A + : R^* 9 $(x) 

Using lattice Fourier transformation with periodic boundary condi- 
tions, i.e., 

I>(p) = (27r)-2 e-'P"$(x)A^ 
$(x) = (27r)-2 ^ e^P^'l>(p)r/^ 

this operator has the following simple form in terms of the lattice 
Fourier transform $(p) of $(x) : 

•^(p) ^ + ^iv) 

vM=0 / 
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Therefore the kernel K{x,y) (the matrix) of the hnear operator —A + 
is: 



pgf4 \n=o 



and the kernel of its inverse is accordingly 

K-\x, y) = i2nr E ^"^^"^^ (tl ^^^^^ + j v'- 

(3.3) 



A2 



Note that (13.31) can be written as 
where 

3 



./ n9 9 v^2 — 2cos»uA 



A2 

Accordingly we calculate and estimate, for x G F and some B ^ 
which later will be chosen to equal A{p), the one dimensional lattice 
sum 

^ (2-2cospA)/A2 + 52''- ^^-^^ 
per 

The result is: 

Proposition 3.1. Assume B ^ and \ argi?| < 7r/4. Then one has, 
for all X 

^ (2-2cospA)/A2 + fi2^- ^^'^^ 
per 



_ 27r(l + Afi[V4 + A252/2 + A5/2])-l^l/^ 

~ SV4 + A252 ' ^ 

with z± given in l{3.9\) . 

If M, N & *N are infinitely large numbers and, in the case B is 
infinitely large, 6 = AB[\/A + A'^B'^/2 + A5/2] is infinitesimal, then 
Eq. ( [J. 7] ) can he continued by 

^ -B[V4+A2B2/2+Afi]|x| 

= , (3.8) 

B^JA + A252 ^ ' 

for some ~ e, which is near 27re^^'^'/2i?. 
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Proof. In order to evaluate the sum fl3.5p we first rewrite it as 

^ (2-2cosj9A)/A2 + 52^ ~ ^ (2-e*P^-e-*f^)/A2 + 52^ ~ 
per per 

(2eipA _ gi2pA _ i)M2 + eip^52^ = (2;2 - _ + ^2^2^^- 
per per 



for z = e*^^. For the decomposition into partial fractions we determine 
the zeros of the denominator, as a function of the complex variable z, 
i.e., of 

z^ - {2 + A'^B^)z + l = 0. 

These zeros are 



2 + A252 ± A5 V4 + A252 ^ ^ 

z = z± = ; (3.9) 



and we can write 

1 



22; — ^2 _ 1 -|- 2A2i?2 Z^ — Z^ \Z — Z- Z — Z. 

Under our assumptions for B we know that Re 2+ > 1, \z_\ < 1 and 
z+ ■ Z- = 1. This allows us to use a geometric series to evaluate the 
lattice sum. Under these conditions we get 



z I — n ^ I — n \ ^ / 



z — Z- z — z+ . 
and accordingly the evaluation of the lattice sum is continued by 

A27^ 

pgf \ fc=0 k=0 



^-fe-lgip(fc+l)A 



z+ — z_ 

^ k=0 



per 




Now observe that lattice points are of the form x = koA for some 
integer ko, —L + 1 < ko < L while points of the dual lattice have the 
form p = jr], —L + I < j < L. For m G Z one has the following cases 



L 



. \2L m = 

per j=-i+i i-e""T ' 
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Accordingly fl3.1Up equals 

A 2 °° 
^ k=0 

z+ — z_ ^ z+ — z^ 

By inserting the expression (13. 9p . the value (13.71) for the one dimen- 
sional lattice sum follows. 

Now assume that M, G *N are infinitely large numbers. Denote 
u = 5/\5\, where 5 = AB[^/IT^B^/2 + AB/2]. Then we have 

(1 + Afi[V4 + A2fi72 + AB/2])-l^'l/^ _ 

[(1 + w|5|)l/|5|]-|5|l-l/A g-B[vTrA^/2+AB/2]|x| 

5 V4 + A252 ~ + A252 

where we put 

(l + n|<5|)i/l^l = e:, 

i.e., 6=1, = 6^'°^^^+"''^'^ and ^log(l ^ 1, if 5 is infinitesimally 

small. □ 

Remark 3.2. Among other things our calculations for the lattice sum 
have established that for x G F, 

per 



The continuum version of this result reads 



e^P^ - Z^ - Z- Z^ I 



dp = = 2ttz 



x/A 



It is interesting to note that the summation and the integration give 
precisely the same value. 

Proposition 3.3. Let M, N be infinitely large numbers and Mq = 
VM. If \p\ < Mo, then for xq G T, 
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and e*=i.(p) ~ e. // |p| > Mq then 



3 ^ -1 



A2 



< (27r)-='2-2M,|:>=olA^^. (3.12) 

Proof. Recall the definition of A{p) in (13 .4^ . A basic estimate for the 
cosine, i.e., < 1 — cost < ^t^ for \t\ < tt, yields 

4 



m'' + —p^ < A{py <m^+p 



2 



since |p^A| < vr for p G F"'. 

We prepare the application of Proposition 13.11 with B = A = A{p) 
by checking that 6 = AA[V4 + AM2/2 + AA/2] is infinitesimal. On 
the basis of the above estimate for A{p)^ this is straightforward for 
\p\ < Mq. Thus (EH]) follows. 

Since Ay/I+^A^ and z+ = {2 + A^A^ + AAV^ + A^A^)/2 are 
mcreasmg functions of A > 0, z'^"""^ / A^IT'^A^ is a decreasing func- 
tion of A and thus is estimated from above by its value at the minimum 
value Aq for A{p) for |p| > Mq = \[M. By our estimate for Aipf' it fol- 



lows Aq < 2\fM and again AAo and 5o = AAq\^ ^ + AM2/2 + AAo/2] 
are infinitesimal. Hence Proposition 13.11 applies and (I3.12p follows from 
the lower bound Aq>\Mq. □ 

Next we prepare the evaluation of the 4-dimensional lattice sum by 
two lemmas. 

Lemma 3.4. Let M, G *N he infinitely large numbers. If Xq E T is 
not infinitesimal, then 

-1 



2 \ _4 

H- in. 

A2 

pef4 \a'=o 



-V 



Proof. Let Mq = yM and suppose that |xo| is not infinitesimal. Then, 
by Proposition 13.31 

-1 



2 — 2 cosp^A 2 

pefMp|>A/opoef Vm=o 



n~V2Mo\xo\/7T q-2VM\xo\/tt 

^(2-)-' E ^TTfTTTr"' S (2')-"=^' 



2y2Mo/7r 22yM/7r 
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(2vr) 



^(p)-VlgF+m2|xoi 



2v/|qP + m2 



<(27r) 



,2-2|Pll^o|/vi 

4|p|/vr 



for any standard e > 0, there exists a finite Mi > such that 



pef3,Mi<|p|<J\/o 



-V 



< e 



if l^ol is not infinitesimaL This shows that for all e > there exists Mi 
such that 



PGr3,|p|>Afi poef \M=o 
We also have 



2 — 2 cosp^A 2 



-1 



A2 



+ m ?7 



< e. 



V e > 3 Afi 



E 



-i/|p|2+m2|xo| 



pefMp|>Afi 



12 I 



< e. 



Let Ml > be finite. If |p| < Mi, then one has, for some < 6^ 
9(Pp) < 1, 

and 



2v^ 



p + 



2v/|p|2 + m2 



It follows that 



(27r) 



P6f3,|p|<Mi 



2v/k|2 + m2 



-(2 



E 



P6rMp|<A/i 



2v/|p|2 + m2 



and hence 



(2-)-' E 



gipa;g-^|p|2+m2|xo| 



2v/|p|2 + m2 



-7/ 



□ 
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Lemma 3.5. Assume that M, N N are infinitely large numbers. If 
X is finite and \xq\ not infinitesimal, then the following lattice sum is 
infinitesimally close to the expected integral: 



(2-)-^ E - (2vr)-^ / /^^^ dp. 



Proof. For 



2v/|p|2 + m2 



2v^ 



2 _|_ 



2 _|_ j^2 



calculate 



^ipx^^^\p\^+m?\xo\ 



^/W+ 



m 



2^3 



^ipx^^^/\p\^+^\xo\ 



2 _|_ ^2 



and estimate 
9 



f{x,p] 



< e 



v/|p|2 + m2 |p|2 + m2 



Therefore the variation of f{x,p) on — ''7/2, + ^7/2] is smaller 

than 



N 

This shows that 



V f{x,p)r]'^- / 



f{x,p)dp 



< 3 



per^ 



iV 



(|a;o| + \x^\)e~V\P\'+"^'\-o\ 



+3^ 



g-A/ipP+^|a:o| 



dp 



<3^(|xo| + |x^|)|e- 



-m\xo\/\/2 



Since 



/ f{x,p)dp^O, 

J\p\>y^M 



dp 

g-|p||xo|/V2 
g-A/ipP+^|2^ol 

-(ip 0. 
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this lemma is proved. □ 

By combining Lemmas 3.4 and 3.5 we arrive at the main result of 
this section. 

Theorem 3.6. Assume that M, N ^ *N are infinitely large numbers. 
If x,y G is finite and \xo — yol is not infinitesimal, then the lattice 
sum ^3. 3|) is infinitesimally close to the expected integral 

K-\x, y) = (2vr)- e-(-^) (j^ ' " '^"^^^^ 



per* \M=o 



2 1 4 

m rj 



-(2vr)-^/ dp (3.13) 

According to this result, the continuum limit of K~^{x, y) is the two 
point Schwinger function Sm{x — y) of the free neutral scalar field of 
mass m. In fact, if x and y are standard real number and xq ^ yo then 



^ip(x-y) ^~^/\p\^+m^\xo~yo\ 

'Smix -y)= I ^ /I I, , , dp 



/• ^ip{x-y)^-y/\p\^+m2\xo-yo\ 
= / o /I 12^ 2 dp = S„,{x-y). 

This follows from the transfer principle of nonstandard analysis and 
means that integrations in both standard universe and nonstandard 
universe coincide. For finite x, ?/ G F such that Xq — is not infinites- 
imal, we have 

Smix-y) = K~^{x,y) ^ *S^{x-y) ^ *5'„(st x-st y) = S'^(st x-st y), 

where stx is the standard part of x, i.e., the unique standard real 
number infinitesimally close to x. The two point Wightman function 
thus is 

lim S^{i{xo -yo) + e,x- y) 



€^ + 



lim Hxo -yo-ie,x-y) = dI^\xq -yo,x- y) 



I- pip{X-y) ^-i^J\p\'^+m?{xQ-yo-ie) 

= hm (27r)-3 / dp. (3.14) 

4. Convergence of the lattice approximation of the two 
POINT Schwinger function for the free Dirac field 

We denote \l/(x) = (\l/i(x), . . . , \l/4(x))^, and recall the notation in- 
troduced in Section 2. The discrete version of the Dirac operator then 
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IS 



E 

/i=0 



+ m. 



(4.1) 



The kernel of its inverse will be the lattice form of the 2-point Schwinger 
function for the free Dirac field (compare (12. 6p ). Naturally, we deter- 
mine this inverse in analogy to the continuum case and use lattice 
Fourier transformation instead of the standard Fourier transformation. 
The lattice Fourier transformation transforms 7^V^\E'(x) into 

(e*p.^_l)/A 

{l-e-'P^^)/A 











if J 



1, 2, 3, with Qj 


-CTq 







^ ajqj 
- ^3(13 
i{l — e~*^^'^)/A, respectively into 

e^P«^-l)/A 

(l-e~^P»^)/A 



#(p) 





0. Thus the Dirac operator 

3 







l'(p). 



is transformed into 



iqo + m cr ■ q 
— (7 q —iqo + rh 

In order to calculate the inverse of the Dirac operator, in analogy to the 
continuum case, we calculate first the inverse of the second order dif- 
ferential operator of which the Dirac operator is a factor. Accordingly 
we determine first this second order operator. 

Under lattice Fourier transformation the operator 

3 



where V' = Pi V: 



P_Vt, is transformed into 



-iqo + m —cr ■ q 
a q iqo + rh 



In order to calculate the composition 



tqo + m 
-crq 



(yq 

-iqo + m 



-tqo + m —cr ■ q 
cr q iqo + fh 



(4.2) 
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of the above two operators we introduce the abbreviations a = iq^ + m 
and b = cr ■ q and find 

f a b \ f a -b \ _ f aa + bb Q _ \ _ [K 0\ 
\-b a J \ b a J ~ \ aa + bb J ~ \ K* ) 

with 

K = aa + bb = (|goP + \q\'^ +Tn^ + 2im{qo - qo))(7o 

+ {qiQ2 - q2qi)<yz + (^293 - qzq2)(yi + {qzQi - qiqz)(y2- 

We decompose the matrix K into its Hermitian and anti-Hermitian 
part 

K = D + 2iE, K* = D- 2iE 

where 

D = (IgoT + + m^ + 2zm(go - qo))(To, 
E = (Im 5253)0^1 + (Img3gi)a-2 + (Im 51^2)0-3. 

Observe that 

Kqo - Qo) = 2Imgo = 2(1 - cospoA)/A > 0, 

qjqk - qkqj = 2ilm {qjqk)- 

Next we calculate the eigenvalues of the Hermitian matrix E from the 
equation 

det(£' - A(To) = - (Img2g3)^ - {^'^q^qif - (Imgig2)^ = 
and find 

A = ±v^(Img2g3)^ + (Imgagi)^ + (Imgig2)^ = ±P- 

There exist orthonormal eigenvectors x± of E, i.e., Ex± = +px±, which 
are also eigenvectors of D and thus of K: 

Dx± = (IgoT + + + 2im{qo - qo))x± = kx± 

and 

Kx± = {D + 2iE)x± = {k± 2ip)x±. 
It follows, for any a± G C, that K{a+x+ + a^x^) = a+{K + 2ip)x+ + 
a-{K, — 2ip)x- and therefore H-R'xll^ = {k^ + 4p^) for all vectors 
X. Therefore {k^ + Ap^)~2K is a unitary matrix. It also follows that 
is a unitary matrix too and thus the following relations 

hold: 

^JH.^ + Ap^K*-' = Kj^^T+li?, K*^' = K/{k^ + Ap^). 

Now it is straightforward to calculate the inverse of the product 
operator fl4.2p : 

f —'iqo + ^ — cr ■ q \ ^ Z' iqo + m cr q 
\ (T q iqo +rh J y - a q -iqo + rh 



— 1 
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K-^ \ _ 1 f K* 
K*-' J ~ K^ + Ap^\ K 

from which the inverse of the Dirac operator in lattice Fourier trans- 
formed form is easily calculated as 



tqo + m cr ■ q 

— a q —iqo + rh 

1 f —iQo + — c ■ Q \ ( K* 



+ Ap^ \ fT q iqo + m J \ K 



(4.3) 



In the rest of this section we are going to calculate the lattice Fourier 
transform of this identity and will show that the continuum limit of 

<-)-E^(A-:-)(Tl-)"' '^-^^ 

equals the well known integral representation of the two-point function. 
As in the scalar case, the four dimensional lattice sum is evaluated 
successively, beginning with the sum Ylipaet- After some preparations, 
a succession of lemmas will prepare the final result of this section. 
Theorem 14.111 

In order to calculate the matrix 



+ Ap'Y^K = {k^ + 4p2)-i(«:ao + 2iE) 



we determine first the factor 



K 1/1 1 

+ 



K2 + 4p2 2\K + 2ip K-2ip 

and expand k = |goP + + + 2im(go — Qo) = (1 — 2mA)(2 
2 cospoA) /A^ +\q\'^ + w?. 

Next we prepare the evaluation of the sum 



E 



2mA ^ (2 - 2 cospA)/A2 + Bl 

per 



5i = jl!±^, A^ = \q\^ + m^ 
^ l-2mA' 



with the help of Proposition 13.11 by the following lemma. 

Lemma 4.1. For the quantities introduced above these statements hold: 
a)0< p< \q\^l\f2, < ±arg5± < 7r/6; 

h) ij \p\ < Mo = VM, then p < 0r|q|7(V2Mo) and arg5± ^ z/M 
is an infinitely large number. 
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Proof. Note that 

qkQjA^ = {e'"^^ - l){e-''"'^ - 1) = e^^f^-^^)^ - e'^^^ - e-'^"^ + 1, 
Im qkQjA^ = sm{pj — p^) A + sinpjA + sinpfcA 
= sinpjA cospfcA — cospjA sinpfcA — sinpjA + sinpfcA 
= sinpfcA(l — cospjA) — sinj>jA(l — cos j»fcA), and 
(Img^gjA^)^ = sin^pfcA(l — cospjA)^ 

— 2 smpfcA(l — cospjA) sinpjA(l — cospfcA) +sin^pjA(l — cosp^A)^, 

pW = 

[(Img2g3)^ + (Imgs^Zi)^ + (Imgig2)^]A'^ = sin^p2A(l - cospsA)^ 

— 2 sinp2A(l — C0SP3A) sinp3A(l — C0SP2A) + sin^ j93A(l — cos j>2A)^ 
+ sin^p3A(l — cospiA)^ — 2 sm]?3A(l — cospiA) sinpiA(l — C0SP3A) 
+ sin^piA(l — cosp3A)^ + sin^piA(l — cosp2A)^ 

— 2 sinpiA(l — C0SP2A) sinp2A(l — cospiA) + sin^p2A(l — cos j>iA)^ 
= (1 — cosj>iA)^[sin^p3A + sin^p2A] + (1 — cosp2A)^[sin^p3A+ 

+ sin^piA] + (1 — cosp3A)^[sin^p2A + sin2piA] 

— 2 sinp2A sinp3A(l — cosp3A)(l — C0SP2A) 

— 2 sinp3A sinpiA(l — cospiA)(l — C0SP3A) 

— 2 sinpiA sinj92A(l — cosp2A)(l — cospiA). 

Since | sinp^A| < 1, we have 

3 

p^A^ < 2^(1 - cospjA)^ + 2(1 - cosp3A)(l - C0SP2A) 
i=i 

+ 2(1 - cospiA)(l - C0SP3A) + 2(1 - cosp2A)(l - cospiA) 

Thus we have p < \q\^ Since < ±argS|., 

< ±arg5^ < tan"^ V2 < 7r/3, < ±arg5± < 7r/6. 

If IpI < Mo, then | sinp^A| < |p^A| < .JiiM^^ , and therefore p^A^ < 
7r(|q|2)2A7(2M2) and p < 0F|g| V( V2Mo). Hence tan(arg5^) < 
-^\PkJ2 and thus argi?-|- ^ if M is infinitely large. □ 

By Lemma 14.11 B± satisfies the conditions for the constant B in 
Proposition 13.11 so this proposition applies for the present case and 
yields 
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Proposition 4.2. For the quantities B± introduced above we have 
B± 7^ 0, I argi?-|-| < 7r/6 and, for all x G F, 



E 



^ (2 -2cospA)/A2 + 5i 
per 



V 



27r(l + A5±[v/4 + A2fi2/2 + A5±/2])-l^l/^ 
fi±V4 + A2fi2= • 



//M, N E *N are infinitely large numbers and 6 = AB±[^/A^^~^^^/2- 
AB±/2] is infinitesimal, then, for some e^, ~ e, the above sum equals 



27re* 



B±[^4+A2b2/2+AB±]|x| 



5±v/4 + A252 

w;/izc/i near 27re~^± 1^1 /25± and /ess i/ian 27r|2-^l^l/25±| . 

Our sum is evaluated further and estimated in the next proposition. 

Proposition 4.3. Let M, N E*N be infinitely large numbers and Mq = 
VM. If \p\ < Mq, then for Xq G T 

ixp 



-V 



:(2vr) 



(2-2cospoA)/A2 + 5| 

ipx^ („]~B±ly/4+A^Bl/2+AB±/2]\xo\ 



poer 



B±^A + A^B. 



(4.5) 



and e^{p) ~ e. // \p\ > Mq then 

(2")-' E 



^ixp 



poGF 

<(27r)^^2~^^'^°''^°''''^ 



(2-2cos|)A)/A2 + fi2^ 



AMo/tt' 

Proof If \p\ < Mo, then p < 0F|g|2/(V2Mo) and 



(4.6) 



^|g|VMo ^ ^2l±VW^ < 



M2 



— 1 + V27r/M, 



1 + A - 1 + A -V'J i + A 

This shows that A5± and 6 = A5±[^4 + A25|/2 + Afi±/2] are 
infinitesimal. Application of Proposition 14.21 implies 04.51) . 

Apply Lemma l4.ll to argi?-|- = 6±. This gives |6'-|-| < tt/G. For 
<P± = arg v/4 + A25| we get < ±0± < ±9± and. 



> 



^2 



2mA 



, |4 + A^^il > 4 + A^ 



since 



A' 



2mA' 
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it follows 



1 - 2m A ' 



Re5±y'4 + A25| = |S±||y^4 + A2B||Ree'^±e''^± 

> -|5±|K/4 + A25i| > i — AV4 + A2A2. 

- 2' ^"V ^ ±1 - 2(l-2mA) 

Since (recall that z± is defined by Eq. (13. 9p with 5 replaced by B±) 



, , 2 + Re A25i + Re A5± a/4 + A25i 
Iz+I > Re2;+ = 



> 1 + (1/2)[A2a2 + (1/2)AAV4 + A2A2]/(1 - 2mA), 

we have 

|^-|^o|/A| 

|5±v/4 + A252| - 



{1 + (l/4)Ay4[V4 + A2A2 + 2AA]/(1 - 2mA)}-l^'ol/^ 

+ A2A2/(l-2mA) ■ ^ ■ ' 



If IpI = Mo, then b = (1/4)AA[V4 + AM2 + 2AA]/(1 - 2mA) is 
infinitesimal, and the right hand side of Eq. (14. 7p is: 



j(l + ^y/5y5\xQ\/A g-{(l/4)^[V4+A^A^+2AA]/(l-2mA)}|xo| 

AV4 + A2A2 ~ AV4 + A2A2 

/ l/(l-2rnA)x_{(i/2)A}|xo| 2-^(V2)A}|xo| 



2A - 2A - 4Mo/7r' 

where we used again the fact that 

Since the right hand side of Eq. (14. 7p is a decreasing function of A, we 
estimate, for |p| > Mq, as follows: 

l^-l^ol/A, . 
1^-^ I < 2-2^-fo|2.'o|/7r. 



|fi±V^4 + A25|| ~ 4Mo/7r' 

This proves (14. 6p . □ 

After these preparations, in a sequence of lemmas, the evaluation of 
the four dimensional lattice sum for the Dirac field is done successively. 
The first step is: 
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Lemma 4.4. For infinitely large numbers M, N G *N, if Xq & T is 

not infinitesimal, the four dimensional lattice sum is approximated by 
a three dimensional one: 



;2-2cosj)oA)/A2 + 5| 



2v/|p|2 + m2 



Proof. With the abbreviation Mq = a/M we estimate as follows: 



(2 



:7r 



pef3jp|>Mopo6r 

2-2A/o|a:o|/7r 



(2-2cosj9oA)/A2 + 5| 



pefMp|>Mo 



4Mo/7r 



2-2A/o|xo|/7r 

< (27r)-^/2jvf3____ ^ 0. 



4Mo/7r 



Since 



5±V4TA25| 

2-A|xo| 2"2|p||xo|/7r 

- 2y4 - A\p\/ti 
for IpI < Mq, for any standard e > 0, there exists a finite Mi > such 
that 



pef^A/i<|p|<A/o 



gipa; - B± [^4+ /2+ AB± /2] I | 



< e 



if IxqI is not infinitesimal. This shows that 



V e > 3 Ml (27r)-' 

We also know 

V e > 3 Ml 



E E 

pef3,A/i<|p| poef 



(2 -2cospoA)/A2 + 



< e. 



(2 



E 



-^|p|2+m2|xo| 



pet''\\p\>Mi 

If IpI < Ml for a finite Mi > 0, then it follows that 



< e. 



and 



p B±K. \Jp'^ + im? 

e* (P) [ V^+^'-B2 /2+ AB±/2] |xo I g-VlPP+™'l^o| 



2v^^ 
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Hence, for non- infinitesimal |xo|, we conclude 



^ ^ ^ipx^^(^pyB±[^i+A^Bl/2+AB±/2\\xo\ ^ 



P6rMp|<A/i ± V ' ± 



and the proof is complete. □ 



Lemma 4.5. Under the same condition as Lemma\4^ the four dimen- 
sional lattice sum 

H 

(2-2cos|)oA)/A2 + 5|' 
is infinitesimally close to the spatial three dimensional one 



(9 ^ gfc 4 



^3 2v/|p|2 + m2 

/or A; = 1,2,3. 

Proof. We can proceed in the same way as in the proof of Lemma 



Lemma 4.6. Under the same condition as Lemma \4-4\ the four di- 
mensional lattice sum for the Dirac operator is reduced to a three di- 
mensional one in the following way: 



4 



E ipx -y/\p\^+fh?\xo\ / ~ rr .■n\ 

^^-3 2v/|p|2 + m2 V '^■P ^ ; 
Proo/. Note that if p 7^ 

^ixopo ^ gixopo 



2 U- 



1 



2ip K + 2ip 



and 



- V e"""'— ( i i ) Er7, 

^ 2p\K-2ip K + 2ipJ " 

Poef 

+ Ap^ ^ K 
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if p = 0. Since 

El 1 . . pixoPo 

K±2ip' 1 + A ^ (2 - 2 cospoA)/A2 + Bl ' 
poef poef 

2tx (1 + A5±[v/4TA25|/2 + A5±/2])I^oI/^ 
~ 1 + A 5±v/4 + A25f 

and E/p is finite, this lemma follows from Lemmas 14.41 and 14.51 □ 



Lemma 4.7. Under the same condition as Lemma \4-4\ the following 
approximation holds: 



^ ' (2 - 2cospoA)/A2 + 52 



pipx ~^/\p\^+^\xO\ 



1x0 



V 



Proof. 

^ (2-2cosj9A)/A2 + 5^^ " ^ (2-2cosj9A)/A2 + 5| 
poef poer 

_ 27rV+^"°/^ _ 27r^!"/^(z_ - 1)/A 

2;+ — z„ z+ — 

_ 27rz!°/^ A2E2 _ ^^^4 + A2g2 

2;+ — 2_ 2 
if > 0, and if < 0, then 

_ 27r4°/^ A2^2 ^ ^ ^2^2 
2;+ — 2;_ 2 
Now we can prove this lemma in the same way as Lemma 14. 4[ since 

A252 ± ^4 + A25| _ 1 
2V4 + A252 2 
for IpI < Mq. □ 



Lemma 4.8. Under the same condition as Lemma \4.J\ the following 
approximation 

^^^^ 2^ ^2 + 4 2 1 % J V ^ 



of the four dimensional by a three dimensional lattice sum holds. 
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Proof. We can prove this lemma in the same way as Lemma 14.51 □ 
The combination of lemmas 14.61 and 14.81 gives 



Lemma 4.9. Under the same condition as Lemma \4-4\ ^^e four di- 
mensional lattice sum for the Dirac operator is reduced to a three di- 
mensional one: 



^ ^ + 4 2 I q + m M K ) 



4 



7] 



Lemma 4.10. Assume that M,N G *N are infinitely large numbers. 
If X is finite and \xq\ not infinitesimal, then the 'spatial' lattice sum is 
approximated by the expected integral: 



Epipx -y/\p\^+m'^\xo\ / rr . n \ 



pgp3 V 1^1 



Proof. The proof strategy of Lemma 3.5 applies. □ 
Using the formulae [A > 0) 

we get the main result of this section. 
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Theorem 4.11. For infinitely large numbers M,N G *N, for all finite 
X G F"^ for which |xo| is not infinitesimal, the lattice Fourier transform 
of Iji4-^ is given by 



iqo + m -cr ■ q \ f K* 



- (2vr)- / ( + ^ .-"^ ^ ) dp. (4.8) 

5. Convergence of the lattice approximation for the 
interacting theory - in the sense of 
ultrahyperfunctions 

For a motivation of relativistic quantum field theory in terms of 
tempered uhra-hyperfunctions as the appropriate framework for a rel- 
ativistic quantum field theory with a fundamental length and for a 
brief introduction to the mathematics of such a theory we have to refer 
to p]. Here we just mention the basic definitions and results about 
tempered ultra-hyperfunctions as we need them. 

For a subset A of M", we denote by T{A) = W +iA C C" the tubular 
set with base A. For a convex compact set K of M", %{T{K)) is, by 
definition, the space of all continuous functions / on T{K) which are 
holomorphic in the interior of T{K) and satisfy 

\\ff(K), = snp{\zPf{z)\; z G T(K), |p| < j} < oo, j = 0, 1, . . . 

where p = (pi, . . . ,Pn) and = zF^ ■ ■ ■ z^ . Tb{T{K)) is a Frechet space 
with the semi- norms H/H"^'-^'''-''. If Ki C K2 are two compact convex 
sets, we have the canonical injections: 

%{T{K,))^%{T{K,)). 

Let O be a convex open set in M". We define 

T{T{0))=\im%{T{K,)), 

where Ki runs through the convex compact sets contained in O and 
the projective limit is taken following the restriction mappings. 

Definition 5.1. A tempered ultra- hyperf unction is by definition a con- 
tinuous linear functional onT{T{W^)). 

Characterizations of tempered ultra-hyperfunctions are known since 
many years ( [HI HEl [IT] ) . The most convenient one for our purposes is 
based on a result in [2] which we prepare briefly. 
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Let ^o(W^) be the space of all functions F, holomorphic in an open 
set W C C"", with the property that for any positive numbers e, K, 
there exist a multi-index p and a constant C > such that 

\F{z)\ < C(l + \zP\) for all z G C"\(C"\iy)„ \lmzj\ < K 

where {'C^\W)^ is the open e— neighbourhood of {<C^\W). Let a = 
((Ti, . . . , (T„) be a vector with components ctj G {±1}. For such a vector 
a and a number k > introduce the open set 

= {(^1) ■■■,Zn) E C"; (Tjim 2j > /c, for j = 1, . . . , n} 

and the space ^o(C^^) introduced above. Next we consider collections 
{Fo-} of elements F^- G ^o(C";.). Furthermore, for e > 0, > 0, and 
G {±1}, define the path 

= ^aji^i /c) == {2 G C; 2 = X + iaj{k + e),x G M} . 

and then the product path To- = YYj=i ■ 

Then from the definition of the spaces T(T(M")) and Ao{C^f^) it is 
clear that for any collection {F„} of Fg- G ^o(C";.) the assignment 

r(T(M"))9/^({F4,/) = Vai---a„, / F.(2)/(2)dz G C (5.1) 

is well defined and for fixed collection {F^} is linear and continuous in 
/ G r(T(M'^)). Thus for given collection {F^}, F^ G 

nT{K))3f^{{F^}J) 

is a tempered ultra-hyperfunction. Conversely, it is shown in [TTl [2] 
that for any element M of T(T(M"))', there exist constant k > and 
a collection {F„} of functions F^- in Ao{C^i.) such that 

M(/) = ({FJ,/) (5.2) 

for all / G r(T(M")) (see also [IHIISlin]). This proves 

Theorem 5.2 (characterization tempered ultra-hyperfunctions). ^ lin- 
ear functional M on T(T(]R")) a tempered ultra-hyperfunction if, and 
only if, it is of the form ( 15. ( 15.^) for some k > and some collection 
{Fo-} of functions F^ in ^o(C" fc)- 

Remark 5.3. In quantum field theory with a fundamental length, often 
functionals appear which are defined for g G T(T(]R^'^)) by 

{F,g) = J F{x^—X2—i{k-\-e),Xi—X2)g{x^—i{k-\-e),Xi,X2,X2)dx^---dxl 

= j F(x° - i{k + e), x)f{x^ - i{k + e), x)dx° ■ ■ ■ dx^ (5.3) 

for an analytic function F G v4o(Vr) defined in the region 

W = {z = {z^,...,z^) G C^-Im2° > \lmz\ + k} 
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and 



j g{z + X2, X2)dxl ■■■dxle T(T(M^)). 



It is clear that the integral (15. 3p defines a tempered ultra-hyperfunction. 
But the integral representation (15.31) looks quite different from the in- 
tegral representation (15.11) . (15. 2p which characterizes tempered ultra- 
hyperfunctions M according to Theorem 15. 2[ Here we explain that 
the integral representation (15. ip . (15.20 can be expressed by the integral 
(15. 3p in certain situations (e.g., the support of the Fourier transforma- 
tion M of M is contained in the forward light-cone Vj^ = {x & Mf^; x^ > 
\x\}). For simplicity, we assume n = 2. Consider the situation that 
= = 0. Then 

(F,/) = 5^aia2 /" F,{z)f{z)dz= Yl ^1^2 /" F,{z)f{z)dz 

cr -^^^ ^=(1 _1)^(_1 _1) -^^T 

F(i,^i){xi+i{k+e),X2-i{k+e))f{xi+i{k+e), X2-i{k+e))dxidx2 

J J F(^~i~i)ixi-i{k+e),X2-i{k+e))f{xi-i{k+e),X2-i{k+e))dxidx2. 
Now we further assume that -F(±i,--i) is analytically continued from 

*^(±i-i),fc = {(^15^2) G C^; ±Imzi > k, -lmz2 > k} 
to the set 

{(^1,2:2) € C^;±lmzi — Imz2 > 2k, — Im2;2 > k}. 
Then by deforming the path of integration, we get 






F{i-i){xi, X2 - i{2k + e))f{xi,X2 - i{2k + e))dxidx2 

J J 

and 





F(„i,i)(xi, X2 - i{2k + e))/(xi, X2 - i{2k + e))dxidx2. 

Put G{zi,Z2) = „i)(2;i, 2:2) - -^(-1,1) (21, 22)- Then G{zi,Z2) is an- 
alytic in 

{(2:1, Z2) e C^; (Im2;i - Im,Z2 > 2k) A (-Im Zi - Im2:2 > 2k)} 
= {(2:1,2:2) e C^; -Im2;2 > |Im2;i| + 2k}, 

and we have 

{F,f)= / G{xi,X2 + i{2k + e))f{xi,X2 + i{2k + e))dxidx2. 
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As we had seen at the end of Section 2, the n-point Wightman func- 
tion of the field ^l:{x) is 

. . . , x„) = (det Ay'/'W'.Jxu ■■■,Xn), (5.4) 

where A is the matrix determined by fl2.9p . i.e., (aj^k), jy k = 1, ■ ■ ■ ,n 
with 

Qjj = 1, Qj^k = flfej = '^hrjhrJ^D^^^Xj - Xfc), if j < k. 

In [2] the functional characterization of a relativistic quantum field 
theory with a fundamental length has been given in terms of six con- 
ditions (RO) ■ ■ ■ (R5). Now we are going to show that the system fl5.4p 
satisfies condition (RO) which states that this systems consists of sym- 
metric tempered ultra-hyperfunctions. The first part of this condition 
(RO) says that the assignment 

is a continuous linear functional on T(T(M^")), for n = 1, 2, 3, . . .. 

In order to investigate this continuity property, we apply the general 
expansion formula for determinants we get 

det A = ^ sgn (j, k,..., I)aija2,k ■ ■ ■ an,i 
= ai,ia2,2 ■ ■ ■ an,n + ^ sgn(j, /c, . . . ,/)aija2,fc ■ ■ -a^,/. 

(j,fc,.. .,07^(1,2,. ..,n) 

Because of the special values of the entries aj^k according to (12.91) we 
see 

det A = 1 + P„(aj,fe) (5.5) 

where Pn{0'j,k) is the sum of homogeneous polynomials of degrees m = 
2, • • ■ , n in the entries aj^k, I < j < k < n with integer coefficients. 
The integral representation for Dm "* as given at the end of Section 3 
easily implies, for every e > 0, the global estimate 

\D^-\x^ -ie,x)\ < {27re)-^ for all x G (5.6) 

It follows that \P{aj^k)\ < 1 if we choose all — y^, j < k, sufficiently 
large and put Zj = (x° + iyj, Xj). Hence for these Zj, {det A{z))~^^'^ = 

(1 + P{aj^k{zj, Zk)))~^^'^ is a bounded analytic function of the xj in a 
tubular domain and therefore, according to Theorem 15.21 

>V;(zi, . . . , Zn) = (det A{z)r'/^W',Jz,, ...,Zn) 

determines a tempered ultra-hyperfunction by the formula 

mf) = I (det A{z)r'l'Wl^{z,, . . . , z^)f{z)dz, (5.7) 

where Tj = + i(?/°, 0, 0, 0), for all / G r(r(M^")), i.e., the first part 
of Condition (RO) is satisfied. 



35 



We conclude that the sequence of Wightman functions satisfies 
G T(r(R^"))' for n = 1,2,3, ... . The second part of Condition 
(RO), i.e., 

>Vfc;\(/*) = m-Xif), r (^1, • • • , ^n) = /(^n,...,^l), 

where -j/^^* = follows easily from the fact that 

Di;;\z,-Zk) = D^-\zk-z,). 

To conclude this section, we have a closer look at the two-point 
function 

Estimate (5.6) shows that \Al'^Dm{xi— X2 — ie, Xi — X2Y\ < 1 if e > ^ = 
l/{V2n), and [1 — Al'^Dmixi — — ie, Xi — 032)^]^^/^ is analytic with 
respect to xi and X2- Then the functional defined by (5.7) for n = 2 
and 1/2 — yi = e > can distinguish the two events only if their distance 
is greater than e (see [2]). Since e > £ is arbitrary, i is the fundamental 
length of our theory. 

6. Conclusion 

The results of this article provide a solution of the linearized version 
of Heisenberg's fundamental equation, on the level of all the n-point 
functionals of the solution fields. This has been achieved by employ- 
ing path integral methods for quantization. In order to have all the 
path integrals well defined and to evaluate them rigorously, a lattice 
approximation was used and the continuum limit of this approxima- 
tion was controlled by using non-standard analysis. This continuum 
limit exists in the framework of tempered ultrahyperfunctions but not 
in the framework of tempered distributions. In this way in particular 
the convergence of the lattice approximations for a free scalar field, a 
free Dirac field and for the interacting fields of this model has been 
established. 

In the second part we are going to show that the sequence of all 
n-point functionals which we have constructed satisfy all the defining 
conditions of a relativistic quantum field theory with a fundamental 
length, in the sense of [2]. We do so by first extending the theory 
of [2] to include scalar as well as spinor fields and then verifying the 
defining condition. In addition we offer an alternative way to calculate 
all the n-point functionals of the theory by use of Wick power series 
which converge in the sense of tempered ultrahyperfunctions. And its 
is shown that the solution fields (0, ip) of (11. 4p - (11.51) can be express 
of a point-wise product 

ij{x) = Mx) ■ e'''"^^"^' : (6.1) 
where ipQ is the free Dirac field. 
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